Phase diagram of a 1 dimensional spin-orbital model 
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We study a 1 dimensional spin-orbital model using both analytical and numerical methods. Renor- 
malization group calculations are performed in the vicinity of a special integrable point in the phase 
diagram with SU(4) symmetry. These indicate the existence of a gapless phase in an extended region 
of the phase diagram, missed in previous studies. This phase is SU(4) invariant at low energies apart 
from the presence of different velocities for spin and orbital degrees of freedom. The phase transition 
into a gapped dimerized phase is in a generalized Kosterlitz-Thouless universality class. The phase 
diagram of this model is sketched using the density matrix renormalization group technique. 

PACS: 75.10.Jm, ll.10.Hi, 11.25.Hf, 75.40Mg 

I. INTRODUCTION 

Spin-orbital models arise in many kinds of mate- 
rials. They b&ve been derived for C@o materialEl, 
LiNiC>2 samplcsEI, and degenerate chains in Na2Ti2Sb20 
compound! In this paper we study a one dimensional 
SU(2) x SU(2) spin-orbital model with Hamiltonian: 



H 



E 



(x + Si ■ Si+i)(y + Ti ■ Tj+i). 



(1.1) 



where Sf and T°" are 5=1/2 spin operators at site i. 
We note that this model has an additional Z 2 symmetry, 
interchanging spin and orbital degrees of freedom, along 
the line x = y. 

Our main conclusion is the groundstate phase diagram 
of Fig. 1. Phases I, II and III and-IV, have been dis- 
cussed extensively in previous worklttil In phase I both 
spin and orbital degrees of freedom are in fully polar- 
ized ferromagnetic states. In phase II the orbital degrees 
of freedom are in the fully polarized ferromagnetic state 
while the spin degrees of freedom are in the standard an- 
tiferromagnetic groundstate and vice versa in phase III. 
Phase IV is a gapped phase with spontaneous dimeriza- 
tion. Our new results concern phase V. The point on 
the V-IV phase boundary at (x, y)=(l/4,ljtt) has SU(4) 
symmetry and is Bethe ansatz integrableB The SU(4) 
symmetry follows from the fact that this Hamiltonian is 
simply a permutation operator, interchanging states on 
neighboring sites. The low energy theory for this model 
is known to be the SU(4)i Wess-Zumino-Witten (WZW) 
model, with central-charge c=3, equivalent to 3 decou- 
pled free bosons .Elftij We show here that the entire ex- 
tended region V is gapless. (Previous work concluded 
that only the Z 2 symmetric line — 1/4 < x = y < 1/4 was 
gapless.) We show that the renormalization group (RG) 
flows in region V are to the IV-V phase boundary line. 
This represents a line of critical points of a rather unusual 
kind. All critical exponents are unchanged along the line 



but there are two different "spin-wave" velocities one for 
spin and one for orbital degrees of freedom. Their ratio 
varies continuously along this critical line. They are equal 
at (x,y) = (1/4, 1/4). At the tri-critical point where II, 

IV and V merge the orbital velocity goes to while the 
spin velocity stays finite. Such a vanishing velocity is 
a natural precursor of a transition into a ferromagnetic 
phase where the dispersion relation is quadratic, rather 
than linear, at small wave-vectors. 

This behavior can be n ujnderstood using non-Abelian 
bosonization techniques.LrH The 5/7(4) i WZW model is 
equivalent to a product of two independent SU (2) 2 WZW 
models, one for spin and one for orbital angular momen- 
tum. The SU(2) 2 WZW models is itself equivalent to 
a triplet of Majorana fermions, a representation of the 
spin-orbital model used by Azaria et al.B However, the 
low energy components of the spin operators with wave- 
vectors near ±7r/2 cannot be represented locally in terms 
of the Majorana fermions whereas they can be so repre- 
sented by the WZW models, making this representation 
more powerful in general. The amplitudes in front of the 
decoupled spin and orbital terms in the SU(2) 2 x SU{2) 2 
WZW Hamiltonian are proportional to the spin and or- 
bital velocities. These are equal for x=y but can be seen 
to be unequal in general. In order to test the validity of 
this picture, we predict the finite size spectrum in region 

V for both periodic and open boundary conditions. This 
takes the general form: 



Ei — Eq 
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[VsXf +W X, 



(1.2) 



where I is the system length. As usual in conformal 
field theory (CFT), the finite size gaps are proportional 
to scaling dimensions of operators corresponding to the 
states. Under the non-Abelian bosonization, each oper- 
ator can be written as a product of a spin and orbital 
operator with the additive scaling dimensions, x*'°. Of 
course, at the symmetric point, v s — v and we recover 
the usual CFT result. As we move around in region V 



the velocities change but the scaling dimensions x^'° do 
not. Thus, we may say that this entire region is SU(4) 
invariant up to a velocity rescaling. 

This type of critical behavior is, strictly speaking, not 
Lorentz invariant and is, in fact, governed by an exactly 
marginal non-Lorentz invariant operator in the low en- 
ergy effective Hamiltonian. However, this type of break- 
ing of Lorentz invariance is essentially trivial and is fa- 
miliar from Tomonoga-Luttinger liquids where the spin 
and charge velocities are different. 

The transition between regions IV and V, along the Z 2 
symmetric line, is in a recently discovered generalization 
of the Kosterlitz-Thouless universality classa character- 
ized by the correlation length (and inverse gap) diverging 
as: 



£ oc exp[yl(x 



-2/3n 



(1.3) 



We verify our RG conclusions, to some extent, by 
density matrix renormalization group (DMRG) work for 
chains of length up to 60. Such numerical verification is 
rendered very difficult by logarithmic corrections which 
we also derive. 

After this work was essentially completed Ref. ( nOI) ap- 
peared on the xxx archive which also discovered the gap- 
less region V. However, there it was concluded from nu- 
merical work that region V actually consisted of 3 differ- 
ent gapless phases characterized by the spin and or bita l 
quantum numbers of the lowest excitation. Our Eq. ( |1.2| ) 
predicts crossing of the lowest excited state as the ratio 
v s /v varies but makes it clear that such a level crossing 
does not correspond to a phase transition. In general dif- 
ferent finite size levels cross at different ratios of v s /v ; 
these crossing points are of no particular significance and 
region V is just characterized by the line of fixed point 
with continuously varying velocity ratio. 

In the next section we discuss non-abelian bosoniza- 
tion of this model and the RG flows, deducing the phase 
diagram. In section III we discuss the finite-size spec- 
trum with both periodic and open boundary conditions 
and in particular, explain the level crossings observed in 
Ref. ( nq). We also calculate logarithmic corrections. In 
Sec. IV we present our DMRG results, corroborating our 
analytical predictions. Sec. V contains conclusions. 



II. NON-ABELIAN BOSONIZATION AND RG 

ANALYSIS 

A convenient way to bosonize this modelUta is to be- 
gin with a generalized Hubbard model at a commensu- 
rate filling where charge excitations are gapped. Upon 
bosonizing the fermion model, it is found that the Hub- 
bard interaction only has the effect of gapping the charge 
bosons leaving various gapless spin-orbital degrees of 
freedom governed by an effective Hamiltonian which is 
conformally invariant up to marginal operators. We be- 
gin by considering the SU(4) invariant case. 



Consider the tight-binding model: 

H = E[(- te I aQc - J+ i + h.c.) + U{c] aa c 



I) 2 



(2.1) 



where c JQa is an electron annihilation operator, j labels 
sites and the repeated spin, and orbital indices, rep- 
resented by Greek and Latin indices respectively, are 
summed from I to 2. We consider the case of 1/4 fill- 
ing, i.e. 1 particle per site. In the large U limit only 
states with exactly one particle on every site have low 
energy, giving the SU(4) invariant version of the spin- 
orbital model with an exchange interaction of 0(t 2 /U). 
The spin and orbital operators are represented by: 



+ a? 

^i Cj H/3a: ■*-% 



c i ~ c iab- 



(2.2) 



The SU(4) symmetric exchange Hamiltonian can be writ- 
ten: 



H = (1/4) 5>/S/ +1 , 



where 



S 



A _ Jaa 



Ci {M )aa c Pb,j 



(2.3) 



(2.4) 



and the Ms are a complete set of 15 4 x 4 traceless 



Hermitean matrices normalized so that: 
trM A M B = {1/2)S AB . 



(2.5) 



The factor of 1/4 was inserted in Eq. (2.3) i n or der 
that the normalization agree with that of Eq. (1.1) at 
x=y=l/4. A convenient choice of these 15 matrices is: 

(^fjt/Vs, 8${a l ) b a /a/8, (Of (o- j )l/V8, (2.6) 

with i,j = 1, 2, 3. Thus the first 6 SU(4) operators 
are the spin and orbital angular momentum operators 
and the additional 9 SU(4) operators combine spin and 
orbital angular momentum. The 15 SU(4) operators, Sf 

are given by the spin and orbital operators S, T and 

We may study the low energy degrees of freedom of 
this model by keeping only Fourier modes of the fermions 
near the Fermi points, ±7r/4. Thus we introduce left and 
right movers, tp, tp: 



\ (Vaa(j>^ /4)J + Vw(j>-'^ /4)j ) . (2.7) 



The hopping term gives, at low energies, a Lorentz in- 
variant free Dirac fermion Hamiltonian density: 



H = iv[^(d/dx)ip - ipi(d/dx)ip], 



(2.8) 



with effective "velocity of light" given by the Fermi ve- 
locity, v = y/2t. 

The Hubbard interaction gives 4 different continuum 
terms upon dropping all oscillating terms. These can be 
conveniently written in terms of left and right charge and 
SU(4) currents: 

J=:V tac Va«:, J A =^ aa {M A f*J ph (2.9) 



where the M matrices are discussed in the prece ding 
paragraph. Using the basis of matrices given in Eq. (2.6) 
we see that the first six SU(4) currents are the spin and 
orbital currents which we write as J\ and 3 % respectively. 
We also define right moving currents J and J A . The 4 
continuum interactions obtained from the lattice Hub- 
bard interaction arc: 



Hint/^v = Xo(J A J A + J A J A ) + Ai (J J + J J) 



-g J A J A +g c JJ. 



(2.10) 



The index A is summed over all 15 values and the cou- 
pling constants are all proportional to U. 

To proceed with an RG asaly-sis of this model it is 
very convenient to bosonizeJIjii3 In order to keep ex- 
plicit track oilhe SU(4) symmetry we use non-abelian 
bosonizationE'El Various fermion bilinears can be repre- 
sented in terms of a free charge boson, (f> c and an SU(A) 
matrix field g. The non-interacting action is a sum of 
the usual free boson action and the WZW action with 
the integer-valued topological coupling constant having 
the value k=l. Both terms in the bosonized free Hamil- 
tonian are quadratic in currents: 

H = (tto/4) (JJ + JJ) + (2irv/5){J A J A + J A J A ). 

(2.11) 

It is a rema rkab le factta that the free fermion Hamilto- 
nian of Eq. (2.8) can also be written in a form quadratic 
in the currents of Eq. ( |2.9[ ). This observation is cru- 
cial for establishing bosonization. We now consider the 
effects of the 4 interaction terms in Eq. ( p. 10 ). The 
non-Lorentz invariant interactions, Ao and Ai renormal- 
ize the amplitudes of the free Hamiltonian corresponding 
to renormalizing the velocities for charge and SU(4) ex- 
citations: 



u 4 ^v(l + 5A ) u c ^v(l + 8Ai). 



(2.12) 



The g c interaction is easily handled since the charge cur- 
rents are linear in the charge boson: 



J = (1/V8ir)d^ c J = (l/V87r)0 + 0c 



(2.13) 



Thus the charge part of the Lagrangian density (in units 
where v c = 1) becomes: 



Cc. 



(l/2)(l-g c /2)(d^ c )(d^ c ). 



(2.14) 



g c can be adsorbed into a rescaling of the charge boson. 
The go interaction is not so trivial but can be seen to 



renormalize to from an initially negative value which it 
obtains for U > 0. Thus, at small U, the low energy the- 
ory for the SU(4) Hubbard model is a type of Tomonoga- 
Luttingcr liquid with decoupled gapless charge and SU(4) 
degrees of freedom. However, as U is increased we expect 
a phase transition into a phase with gapped charge ex- 
citations. This was recently confirmed by T — Monte 
Carlo worko In the continuum description, this transi- 
tion is driven by an Umklapp term which is of 8 order 
in the fermion fields: 

^Umklapp « V' tl2 V' t2 V t2 V tll ^ii^i2^2i^22 + h.c. 

(2.15) 

Such an interaction is generated at order U 2 . It is irrel- 
evant at small U, being of scaling dimension 4. Under 
bosonization it can be expressed as a pure charge oper- 
ator, cos4y / 7r0 c . However, as U is increased, the scaling 
dimension of this operator decreases due to the rescaling 
of 4> c produced by the g\ interaction. It is expected to be- 
come relevant at a critical value of U and produce a gap 
for charge excitations. Importantly, the effective Hamil- 
tonian for the SU(4) degrees of freedom is expected to 
remain the gapless SU(4) WZW model with a marginally 
irrelevant coupling constant go of 0(1). The SU(4) ve- 
locity parameter cannot be determined exactly by this 
bosonization approach but it value is known from the 
Bcthe ansatz solution, v = ir/8. 

While this method of deriving the SU(4)i represen- 
tation of the SU(4) chain is perhaps most familiar to 
condensed matter physicists it can be also be done more 
elegantly by simply projectirig-.out the charge degrees of 
freedom of the free fermions.EH Either way, the result is 
that the low energy degrees of freedom of the SU(4) ex- 
change model can be represented as: 



S- 



(J' 



r 



const[e l( * /2) hr{gM A ) + h.c 



(2.16) 



where g is the fundamental unitary 4x4 matrix field 
of the WZW model. This has scaling dimension 3/4, 
so the correlation function decays with power 3/2. As 
pointed out by Azaria et al.,u another term, oscillating 
at Ak,F — 7r should also be gener ated by higher order 
processes and be included in Eq. ( 2.16 ). In the WZW 
representation, this operator is the dimension 1 primary 
transforming under the (6,6) representation of SU(4)l X 
SU(4)r. (The 6-dimensional representation of SU(4) is 
the 2 index antisymmetric tensor representation; i.e. its 
Young tableau has 1 column and 2 rows.) In our rather 
cumbersome notation, we may write this tensor as: 



$ 



{aa,f3b} 
{jc, Sd} ' 



(2.17) 



where the indices in curly brackets are antisymmetrized. 
The additional term in Eq. (|2.16D then takes the form: 



(-i)wt.$g:2(Ar 



Sd 
1 0b ' 



(2.18) 



In order to study the effects of SU(4) symmetry break- 
ing, it is convenient to use a different but equivalent non- 
abelian bosonization of the SU(4) exchange model. We 
may replace the 5/7(4)i WZW model by a sum of 2 de- 
coupled 51/(2)2 WZW models representing spin and or- 
bital degrees of freedom. The subscript 2 implies that 
the topological coupling constant takes the value k = 2. 
One way of arriving at this result is by using a different 
non-abelian bosonization of the SU(4) Hubbard model in 
which the fermions are represented by the charge boson 
plus these 2 WZW models. Alternatively we may use the 
conformal embedding of this sum of WZW models into 
the 5/7(4)1 model. The validity of this ^presentation 
can be checked from the fermion identity:ta 



J A J A /5 = (JlJl + «/4. 



(2.19) 



Thus the Hamiltonian of the non-interacting model is 
written as a sum of terms quadratic in spin and orbital 
currents only. The SU(4) matrix field g can be replaced 
by product of SU(2) matrix fields representating spin 
and orbital degrees of freedom. These both have scaling 
dimension 3/8 in the SU{2) 2 WZW model so that their 
product has the correct dimension 3/4. Similarly the 
components of $ can be expressed as sums of components 
of the spin-1 primary fields of the SU(2) 2 models. 

Now consider the effect of breaking the 5/7(4) symme- 
try down to 5/7(2) x 5/7(2), corresponding to the general 
spin-orbital model. We use the field theory approach, 
adding small anisotropic exchange terms , ta king the pa- 
rameters x, y in the Hamiltonian of Eq. (1.1) to have the 
values 1/4+Sx and 1/4+Sy. This leads to the anisotropic 
corrections to the continuum limit Hamiltonian: 

STi/v = -[Sx(J s ■ J s + J s - J s ) + Sy{J ■ Jo + Jo- Jo)} 



+const[SxJ s ■ J s + Sy J ■ J ] 



(2.20) 



The constant factor is non-universal but we expect it 
to be positive as can perhaps be seen most simply from 
the Majorana fermion analysis of Azaria et al.H discussed 
below. 

We first consider the effect of the Lorentz invariant 
terms. Including the large go term already present at the 
SU(4) point, the interaction Hamiltonian becomes: 

H [llt /(4Trv) = g J A J A + gJo ■ X + 52 J s • L (2.21) 

Here the bare coupling constant, go has an unknown neg- 
ative value of 0(1) while: 



51 « ASx, g 2 « ASy, 



(2.22) 



for a positive constant, A. 

The corresponding weak coupling (3 functions can be 
derived using standard methods. The simplest method 
is perhaps to just treat these interactions perturbatively 
in the free fermion theory. This is valid because the gap- 
ping of the charge boson by the large Hubbard interac- 
tion doesn't effect the renormalization of the coupling 



constants in the SU(4) sector or the theory. The result 
is: 



dgo 

dl 
dgx 

dl 

dm 

dl 



4<?o+4.go(3i + <?2), 
-2.go.92 + 2#i, 
-2soSi +2g2- 



(2.23) 



To understand the solution of these RG equations, it is 
convenient to first consider the Z 2 symmetric case, g\ = 
<?2, where they reduce to: 



l-rr =^9o+ 8 9o9i, 
l-^- = -2g gi + 2gl 



(2.24) 



Note that along the Z 2 symmetric line in region V of the 
phase diag ram g± = g 2 < 0. In this case, the first of 
Eq. ( p. 24 ) implies that go increases (i.e. its magnitude 
decreases since it is < 0). The second of these equations 
implies that g\ initially decreases (increases in magni- 
tude) since initially \go\ » \gi\- This continues until 
go = g\ after which g\ begins to increase towards as 
does go- It is also instructive to notice Eq. (2.24) imply: 



1 d In g 



= 2(g + 2 9l ) = 



dln(g +gi) d\ng 1 



2 din I 
Thus, along any RG flow: 

\gi(go + gi)\ 2 \go\ 



dlnl 



constant. 



dlnl 



(2.25) 



(2.26) 



This implies, that at long length scales, in region V, 
\go\ oc g\. Thus, even though the bare value of |<7o| is 
much larger than the bare value of \g\\ this situation 
eventually reverses during the RG flow towards cou- 
pling. Some RG trajectories are shown in Fig. 0. Both 
the non-monotonic flow of g\ and the fact that, asymp- 
totically |<7i| >> |<7o| will have important consequences 
for logarithmic corrections to finite size scaling, discussed 
in the next sections. We emphasize the remarkable fact 
that even though the 5/7(4) symmetry is broken down to 
5/7(2) x 5/7(2), this symmetry breaking is irrelevant and 
the full 5/7(4) symmetry still appears in the low energy 
behavior. 

On the other hand for the Sx ~ <?i > case, the RG 
flow runs away from the critical point go = gi =0. This 
implies that the system has an energy gap and a finite 
correlation length £. We can see the universal behavior 
of th e gap gener ation by integrating the RG equation 
( |2.24 ) using Eq.( [2.26 ) to solve for g - Let us consider 

an initial condition go = 0(1) and g± ~ 8x > 0^ This 

implies the constant in the right hand side of Eq. ( 2.26J ) is 
proportional to Sx 2 . The correlation length £ corresponds 
to the scale which makes the running coupling constants 
diverge 



ln£ = 



' dgx — 

Sx 2g\ - {CSxf - CSx^2gf - (CSx) 2 



5x~ 2 '\ 



(2.27) 



where C is a non-universal constant. In this way we 
obtain the behavior of the correlation function or inverse 
gap as a function of Sx: 



A ~ exp[-A(Sx)' 



(2.28) 



corresponding to a type of generalized Kosterlitz- 
Thouless behavior discussed extensively in Ref. ( H). In 
this case, v — 2/3. 

Next, we study the Z 2 asymmetric model with x 7^ y 
and hence g\ ^ g 2 . Noting that the last 2 equations of 
Eq. (2.23) are equivalent to: 



, d(gi+g 2 ) „ , , , 2 . 2^ 

I J, = -2.9o(.9i + 32) + 2(3! + g 2 ), 

,d{$i-92) 



dl 



^9i-92)[(9i+ 92) -9o], (2.29) 



We see that go will still renormalize towards as long as 
initially 31 + 32 < 0. In this case, g\ + g 2 initially flows 
away from but again turns around and flows towards 
when g\ + g 2 = (70(31 + 32)- Similarly, the asymmetry 
parameter, g\ — g 2 initially increases in magnitude but 
eventually also flows to 0. It can be seen that for Sx + 
Sy > 0, the RG equations predict the development of a 
finite gap. Thus the phase boundary between regions IV 
and V should be at at 45° angle to the x-axis, as drawn 
in Fig. I 

From considering only Lorentz invariant operators in 
^int we wou ld include that Z 2 asymmetry is completely 
irrelevant in region V since it doesn't alter the flow to- 
wards the SU(4) symmetric critical point. However, we 
must also cons ider the non-Lorentz invariant interaction 
terms in Eq. (2.20). It is at this point that it becomes 



very convenient to use the SU(2) 2 X SU{2) 2 WZW model, 
which is equivalent to the SU(4)i model as discussed 
above. We then see that the non-Lorentz interaction 
terms just renormalize the coefficients in front of the 2 
terms in the non-interaction Hamiltonian: 



Ho 



[Jo • Jo 



irv s 



(J, 



*J o ' Jo) 



(2.30) 



For small Sx and Sy, the shift in velocities from the SU(4) 
value is linear: 



717 

717 



BSx + CSy 
BSy + CSx. 



(2.31) 



The positive constants B and B are not universal due 
to renormalization of these non-Lorentz invariant terms 
in the effective Hamiltonian by the Lorentz invariant 
terms discussed above. These velocities are, of course, 



equal along the Z 2 symmetric line x — y. Thus, with 
Z 2 symmetry, all breaking of SU(4) symmetry is irrele- 
vant. However, also breaking the Z 2 symmetry produces 
a marginal line of fixed points, which we may regard as 
the phase boundary between regions IV and V in Fig. 1. 
Along this line the critical theory is the SU(2) 2 x SU(2) 2 
WZW models with unequal v s and v . All critical expo- 
nents are constant along this line. The breaking of 517(4) 
symmetry is thus of a very trivial kind. Nonetheless, it 
has important consequences for the finite size spectrum, 
as discussed in the next section. As discussed in the in- 
troduction, the transition from phase IV to phase II or 
III is naturally associated with the vanishing of v or v s 
respectively. 

We remark that the same conclusions can be reached 
by employing the Majorana fermion representation of the 
spin-orbital model used by Azaria et al.H. The 5C/(4)i 
WZW model is equivalent to 3 free bosons or 6 free Ma- 
jorana fermions. The two SU(2) 2 factors are each equiv- 
alent to 3 free Majorana fermions, which transform as 
vectors under the SU(2) symmetry. Breaking the Z 2 
symmetry just gives different velocities to the spin and 
orbital Majorana fermions. The WZW formulation of the 
problem is somewhat more natural than the Majorana 
fermions because all operators in the underlying lattice 
model can be expressed locally in terms of WZW fields. 
This is not true in the Majorana fermion representation. 
The 2kp components of the lattice spin operators are 
non-local in terms of Majorana fermions. This can be 
understood by using another equivalence. Each Majo- 
rana fermion theory can actually be considered to be an 
Ising model. This contains order and disorder opera- 
tors of dimension 1/8 which cannot be expressed locally 
in terms of the Majorana fermions. The SU(2) 2 model 
is equivalent (locally) to a product of 3 Ising models. In 
particular, all components of the fundamental SU(2) ma- 
trix field of the WZW model, of dimension 3/8 can be 
expressed as various products of the 3 Ising order and 
disorder fields. The WZW model is a more natural for- 
mulation of this problem than the product of Ising models 
since it makes the SU(2) symmetries manifest. 

We also remark that the spin-orbital model provides a 
rare example of a lattice model whose critical theory is 
given by WZW models with central charge k > 1. These 
also occur as critical theories for special iategrable SU(2) 
spin chains of higher spin (with k — 2s) H3 However, the 
integrable models represent multi-critical points. Generic 
spin-s Hamiltonians always renormalize to the k=l WZW 
model for half-integer s or develop a gap for integer s. 
The fact that k > 1 WZW models represent unstable 
fixed points can be understood from an RG point of view. 
They contain relevant operators allowed by symmetry in 
the context of spin-s Heisenberg models. However, the 
structure of the spin-orbital model is such as to forbid 
any relevant operators on the 577(2)2 X SU(2) 2 fixed line. 
This extra symmetry can be traced back to the SU(4) in- 
variant model. Translation by 1 site corresponds to mul- 
tiplying the fundamental SU(4) matrix field by a phase 



i7T/2 



Consequently, a produce of 2 of the these fields, 



which gives the antisymmetric tensor primary field of di- 
mension 1, gets multiplied by a minus sign under trans- 
lation. Thus both these operators are forbidden from the 
effective Hamiltonian by translational symmetry. [They 
are not otherwise forbidden since it is possible to make 
diagonal SU(4) singlets from both of them.] Using the 
SU(2)2 x SU(2)2 representation, the antisymmetric ten- 
sor field of the SU(4) model becomes the two symmetric 
tensor fields of the two SU(2) models. These must also 
transform with a minus sign under translation and hence 
are forbidden, unlike the situation for an s=l Heisen- 
berg model where this operator is allowed in the effective 
Hamiltonian and produces a gap. 



III. FINITE SIZE SPECTRUM 

The analysis of the previous section permits a straight- 
forward prediction of the finite size spectrum of this 
model. This is useful for comparing to numerical sim- 
ulations in order to test our conjectured phase diagram. 
The spacing of energy levels vanishes as 1// as the sys- 
tem size I is increased. The coefficients of 1/1 give scal- 
ing dimensions of operators corresponding to the various 
states. Marginally irrelevant operators give corrections 
to the finite size spectrum which vanish as 1/lYa.l. These 
must be taken into account to obtain reasonable agree- 
ment with numerical data for system sizes that are less 
than exponentially large. We first discuss the finite size 
spectrum ignoring marginal operators, obtaining only l/l 
terms. Logarithmic corrections are discussed at the end 
of this section. We only consider the case where the num- 
ber of sites is divisible by 4. The generalization to other 
chain lengths is straightforward but tedious. 

We begin by discussing the SU(4) symmetric model. 
The groundstate energy for a system of length I with 
periodic boundary conditions (PBC) has the form: 



Eq = e l — ttvc/61 



(3.1) 



where the central charge, c = 3. The finite size energy 
levels with periodic boundary conditions are given by: 



Ek -E = (2Trv/l)xi, 



(3.2) 



where Xi is the RG scaling dimension of the operator 
corresponding to the excited state i. The lowest excited 
states for the SU(4) invariant model correspond to the 
fundamental operator, g of the SU(4)i WZW model with 
x = 3/4. This transforms under the (4, 4) representation 
of SU(4) L x SU(i) R . This full chiral SU(4) is broken 
by marginal operators, discussed below. Only the di- 
agonal SU(4) subgroup is an exact symmetry of the lat- 
tice model. Under this subgroup the (4, 4) representation 
decomposes into the adjoint and singlet representations 
(16=15+1). The Hcrmitcan conjugate operator, trans- 
forming as (4, 4) corresponds to another set of states. 



These two sets of states have crystal momenta ±7r/2 re- 
spectively. Thus there are 30 exactly degenerate lowest 
excited states and two higher singlet states which are 
degenerate up to log corrections. These states all have 
crystal momentum ±7r/2. The next lowest energy ex- 
cited states transform under the (6, 6) representation of 
SU(4)l x SU(4)n, corresponding to the primary field <& 
of dimension x=l discussed in Sec. II. Under diagonal 
SU(4) (6,6) decomposes as: 



36= 1 + 15 + 20. 



(3.3) 



Thus we obtain singlet, adjoint and 20 representations all 
with same energies up to log corrections. The 20 repre- 
sentation is real and has a Young tableau with 2 columns 
and 2 rows. These states occur at crystal momentum tt. 
We now consider the effect of breaking the SU(4) sym- 
metries down to SU(2) x SU{2) spin and orbital symme- 
tries, by giving different velocities to the spin and orbital 
parts of the energies of these excited states. We represent 
the quantum numbers of these states by (S, T) the spin 
and orbital angular momentum of the state. The two sec- 
tors both have central charge c=3/2, so the groundstate 
energy is: 

E = e l - (7r/6J)(3/2)(«, + v ). (3.4) 

The adjoint representation decomposes into 

15 -» (1,0) + (0,1) + (1,1). (3.5) 

The x=3/4 5(7(4) i fundamental field is written as a prod- 
uct of x=3/8 fundamental fields of the SU{2) 2 WZW 
models. Thus equal portions of the energies of these 
states come from the spin and orbital sector. Conse- 
quently the energies are all equal for the (1,0), (0,1) and 
(1,1) states at momenta ±7r/2: 



E-E = (2ir/l)(v s +v o )(3/8). 



(3.6) 



The same is true for the singlet state with logarithmically 
higher energy. On the other hand, the 6 representation 
of SU(4) decomposes into (1,0)+(0,1). Thus the various 
components of the (6, 6) primary field become: 

[(S L ,T L ),(S R ,T R )] = [(1,0), (1,0)], [(0,1), (0,1)], 

[(0,l),(l,0)]and[(l,0),(0,l)] (3.7) 

primary fields in the SU(2)2 x 517(2)2 representation. 
The spin 1 primary field of the SU (2)2 WZW model has 
x=l (left and right scaling dimensions 1/2). Thus we see 
that for the 4 different sets of states listed in Eq. (|3.7|), 
the spin and orbital dimensions are: 

(x s ,x o ) = (l,0), (0,1), (1/2,1/2), (1/2,1/2) (3.8) 

respectively. The energies of the various states are given 
by Eq. (1.2). Finally we may decompose these states 
with respect to the exact diagonal SU(2) x SU{2) sym- 
metry. Note in particular that there are states with quan- 
tum numbers (S, T) — (2,0) and (0,2) with energies 



Ei — Eq = 2ttv s /1, and 2ttv /1 



(3.9) 



respectively. 

Close to the Z 2 symmetric case, v s « v , the lowest 
states all have energy (2tt/1)(3/8)(v s + v ) up to log cor- 
rections. These states include the (S,T) — (1,1) state 
which is in fact lowest due to log corrections. On the 
other hand, as the ratio v s /v is decreased a level cross- 
ing eventually occurs and the lowest states have energy 
(2-KVs/l). These states include the (S,T) = (2,0) state 
which is in fact lowest due to log corrections. This occurs 
when 



(3/8)(v s + v ) = v s , 



(3.10) 



i.e. v s = (3/5)w . Thus we expect the quantum numbers 
of the lowest energy state with periodic boundary condi- 
tions to change from (S, T) = (1, 1) to (2, 0) as we move 
away from the Z 2 symmetric line, x = y. Note that this 
level crossing will occur along a line in the (x,y) plane 
which is a finite distance away from the Z 2 symmetric 
line x = y. Exactly such behavior was observed in the 
finite size spectrum with PBC in Ref. ( [!(]). However, 
the present theory makes it clear that this is certainly 
not a phase transition line. All critical exponents have 
the same value on both sides of this line. The same set 
of low energy states occur on both sides. The energies 
are merely shifted by the ratio v s /v . 

We now consider the case of open boundary conditions 
(OBC). Much longer chains can be studied with OBC us- 
ing DMRG and our results presented in the next section 
are all for OBC. In this case the groundstate energy is 
given by: 



E = e l + ei - (7r/240(3/2)(«, + v„). 



(3.11) 



A non-universal boundary energy, ei, appears and the 
l/l term is reduced by a factor of 4. The energy gaps are 
given by: 



Ei -E Q = (ir/l)(v s x l s + v xl). 



(3.12) 



The prefactor is smaller by 2 than in Eq. (1.2) for the 
PBC case. More importantly, the dimensions, x are 
different in this case, corresponding to chiral conformal 
towers taken from the left moving sector only. For the 
SU(A)i WZW model only the conformal tower of the 
identity operator occurs for OBC. Thus the lowest ex- 
cited states transform under the adjoint representation 
of SU(4) and correspond to the 15 chiral current op- 
erators J A . As discussed in the previous section these 
decompose into pure spin or orbital operators and mixed 
operators. The mixed operators have (S,T) = (1, 1) and 
(x s ,x ) = (1/2,1/2) corresponding to a product of the 
(chiral) spin 1 primary operators of the SU(2)2 WZW 
models. Thus we obtain the energies: 



#(1,0) ~ £(o,o) = tv«// 
£(o,i) - ^(0,0) = kv /1 
£(i,i) ~ £(o,o) = (n/2l)(v s + v ), 



(3.13) 



where £(o,o) refers to the groundstate energy. The next 
excited states have x = 2 but we do not consider them 
here. One of the above states will always have lowest 
energy. A level crossing occurs at the Z 2 symmetric line, 
v s = v„. Note that: 



£ 



(i,i) 



[E { 



1,0) 



£ 



(0,1 



]/2 



(3.14) 



ignoring log corrections. We will test this relation nu- 
merically in the next section. 

Finally we consider log corrections, coming from 
marginal operators. Note that these take quite a differ- 
ent form at the SU(4) invariant point, (x,y) — (1/4, 1/4) 
than anywhere else in region V. This is because the only 
marginally irrelevant coupling constant is go defined in 
Eq. ( 2.10| ) in the SU(4) symmetric case. On the other 
hand, when the SU(4) symmetry is broken, even in the 
case where the Z 2 symmetry is preserved, gp fl ows to zero 
much faster than g± ss g 2 defined in Eq. (2.21) according 
to our analysis of the j3 functions in the previous sec- 
tion, as illustrated in Fig. @. Thus, in this case, we may 
asymptotically ignore go and consider only g\ = g 2 which 
gives different log corrections than does go- At interme- 
diate lengths we may expect very complicated finite size 
corrections corresponding to the RG flows discussed in 
the previous section. In particular, we might expect log 
corrections due to go out to some crossover length and 
then log corrections due to g\ — g 2 for longer lengths. 
Furthermore, the fact that \g\\ and \g 2 \ may initially in- 
crease before eventually decreasing, as illustrated in Fig. 
0, means that the amount of breaking of SU(4) symme- 
try may at first appear to increase with increasing length 
before eventually starting to decrease. Evidence for such 
behavior is presented in the next section. 

Log corrections away from the SU(4) point, at long dis- 
tances, follow immediately from earlier work on the ordi- 
nary s=l/2 Hciscnbcrg model since the marginal opera- 
tor, J s ■ J s , is the same one that occurs in that caseliMS 
(Leading log corrections are independent of the Kac- 
Moody central charge, k.) We must merely add the cor- 
rections twice in the spin and orbital sector, with the 
corresponding velocities inserted. 

We first consider PBC. The log corrections to the en- 
ergy gaps for any Virasoro primary states with left and 
right spin Sl and Sr and left and right obital angular 
momentum Tl and Tr are: 



SE 



2tt 
Ttal 



[v s Sl ■ Sr + v T L ■ T R ] 



(3.15) 



In particular, some of the lowest energies, including log 
corrections, are given by: 

2tt, . /3 1 

E (^)- E (m = -(vs+Vo) 8"4ln7 



£(i,o) - £(o,o) 



2tt 

T 

2tt 



1 



£(2,o)-£(o,o)=y^(l lu/ 



8 41nZ 
1 



41nZ 
(3.16) 



Note in particular that, for v s — v , the lowest excited 
state has quantum numbers (1,1). We also see that the 
value of v s /v at which the (2,0) state becomes the lowest 
excited state depends somewhat on /, approaching 3/5 at 
large /. 

We now turn to OBC. In this case, states are classified 
by only a single spin and orbital quantum number: S and 
T. The general log corrections for excited states are given 
by: 



SE 



2/ In/ 



[t>.5(5+l)+t> T(T+l)]. (3.17) 



The energies of the lowest excited states, including log 
corrections are given by: 



-^(1,0) - -^(0,0) 



IT ( 1 



■K I 1 

£(o,i) - £(o,o) = jv I 1 



/ 



In/ 



(3.18) 



Note that for v s = v , (1,1) is the lowest excited state. 

Finally we give log corrections at the SU(4) point, 
(x, y) = (1/4,1/4). For PBC the general formula can 
be written: 



5E = 



I In I 



qA qA 



(3.19) 



a straightforward generalization of the SU(2) case. We 
may write: 

s£si = (i/2)[(s£ + s£)(s£ + si) - sfsf - s£s£] 

= (1/2)[C - C L - C R ], (3.20) 

where C, Cl and Cr are the quadratic Casimir invariants 
for the diagonal SU(4) group and for 5/7(4) L and 5/7 (4) fl 
groups. 

For OBC the corresponding formula is: 



4/ In/ 



S A S' 



4/ In/ 



C. 



(3.21) 



The quadratic Casimir is C = 4 for the adjoint represen- 
tation. Thus the lowest excited states have energy: 



£-£o 



TTV J 1 

T v hi7 



(3.22) 



The velocity, in this case, is determined from the Bcthc 
ansatz solution to have the value, v = n/8. 

We emphasize that these formulas only apply for the 
exactly SU(4) invariant model. For any other points in 
region V we must use the formulas of the previous para- 
graphs. Despite the fact that SU(4) symmetry breaking 
is irrelevant (up to a velocity rescaling) it acts, in a sense, 
like a relevant perturbation as far as log corrections are 
concerned. This simply reflects the fact that the RG flow 
approaches the SU(4) invariant fixed point along a differ- 
ent universal trajectory when SU(4) is broken than when 
it is unbroken. 



IV. DENSITY MATRIX RENORMALIZATION 
GROUP RESULTS 

We now discuss our DMRG calculations Ed We keep 
m = 1100 states in DMRG calculation and the biggest 
truncation error is 5 x 10 ~ 6 . We calculate for chains with 
open boundary conditions with length up to / = 60. 

Our DMRG calculation is a continuation of previous 
studies on this modelJiS'Ei We will use the ground state en- 
ergies, which are most accurately calculated by DMRGcEl, 
to sketch the phase diagram. As pointed out in Sec. II, 
the phase transition between phase IV and V is of in- 
finite order. Thus the singularity of the ground state 
energy e (x, y) as function of parameter (x,y) in Hamil- 
tonian (1.1) is not obvious. The second derivative of 
e(x, y) goes to zero roughly exponentially when we move 
close to the transition point x — y = 1/4 from the dimer- 
ized phase along the x = y line. We use this property of 
the infinite order phase transition to sketch the boundary 
between the gapless phase and the dimerized phase. The 
phase boundaries into the ferromagnetic phases II and 
III are determined by a level crossing. The final phase 
diagram is drawn in Fig. |l|. 

We now discuss other aspects of our results for region 
V. At the SU(4) symmetric point, (x, j/)=(l/4,l/4), the 
lowest excited states are a degenerate set with quantum 
numbers (S,T)=(1,1), (1,0) and (0,1). This energy gap 
is compared to the prediction of Eq. ( p. 22 ) in Fig. 0, 
showing good agreement. 

When we move away from the 5/7(4) point into region 
V along the x — y line of the phase diagram, the energy 
of the (1,0) excitation will have a finite gap above the 
energy for of the (1,1) excitation. Let's call this energy 
difference AE. We will measure AE by the smallest en- 
ergy excitation in the system. We use the energy of the 
(1,1) excitation as the energy unit for AE. AE repre- 
sents roughly the distance of the system from the 517(4) 
invariant model and so is like g\ — g^ in Fig. @. When 
AE = 0, our system is 5/7(4) invariant and g\ = g% = 0. 
When we move farther away from the 5/7(4) point along 
the x = y line in the 5/7(4) gapless phase, AE increase 
and so does g\ = g^- In the RG flow diagram Fig. H, 
we see that when bare g\ = (72 is small for short chain, 
.9i = .92 will first flow away from go axis. So we will 
also observe AE increase when chain length is small for 
(x,y) close to the 5/7(4) point (1/4,1/4). We demon- 
strate these RG properties in Fig. [I|. 

It is difficult to obtain data for long chainsto demon- 
strate the logarithmic scalings show in Eq.( ft.l8 ). Due 
to the complex RG flow as we show in Fig. |2| and Fig. 
[|, very long chains are needed to study this large / scal- 
ing. For the Zi symmetric 5/7(2) x 5/7(2) model (where 
x = y), we have demonstrated the RG flow numerically 
in previous paragraph. 

For the HOT1-Z2 symmetric case, it is even more difficult 
to verify phase V numerically. We show some evidence to 
support that there is a gapless phase in the central part 



V of the phase diagram in Fig. [j] by fitting simple l/l 
and 1/logZ scaling. For an arbitrary point in the SU(4) 
gapless phase x — 1/4 and y = 1/8, we plot (#(i,j) — 
#( ,o))/ sin^) vs. 1/lnZ in Fig. |. The fitting line 

in the figure shows that (#(i,j) — #(o,o)) ~> f f° r our 
calculation indicating that the phase V is gapless. We 
point out that there is no singularity in the ground state 
energy e(x,y) at the symmetric line x = y in the SU(A) 
phase either. In the extended gapless phase V, there are 
two spin velocities v s and v as discussed in the previous 
sections. In Fig. pL we have drawn 



v s (l) = (E {lfi) - £ (0)0 ))/sm(y— -), 

v {l) = (#(o,i) ~ #(o,o))/sin(— j- ), 

y (i,i)(0 = ( E (i,i) - ^(o,o))/sm(y— ). 



We have also drawn 

v s + v (#io 



#oo) + (#oi - #oo) 



8r 



li 



8(#n — #oo) 



(4.1) 



(4.2) 



in Fig. [|. We see that the p redic ted relation #n = 
(#io + #oi)/2 which, from Eq. ( 3.18 ), should be true up 
to log corrections, holds quite accurately for all I shown. 
The numerical results show that the difference between 
#n and (#io + #oi)/2 is less than five percent in our 
short chain calculation. 

In summary, we studied the SU(2) x SU(2) spin-orbital 
model. By a renormalization group study around the 
51/(4) point we showed there is an extended gapless 
phase both in the Z^ symmetric model x — y and the 
asymmetric model x ^ y. This phase has 5'[/(4) sym- 
metry at low energies after rescaling the different spin 
and orbital velocities. The critical theory may also be 
viewed as a product of 517(2)2 WZW models for spin 
and orbital degrees of freedom. We studied the phase 
transitions from the gapless antiferromagnetic phase with 
S#(4) symmetry to the dimer phase and to the ferromag- 
netic phases. 

We used the ground state energy calculated by DMRG 
to sketch out the phase diagram. DMRG calculation also 
shows evidence for an extended gapless phase by giving a 
leading l/l finite size scaling of the gap for chain length 
up to 60 and verifying other aspects of the RG predic- 
tions. 

NOTE ADDED: After, posting this preprint another 
preprint by Azaria et al.E3 appeared that has much over- 
lap with this one. It analyses the complicated scaling 
of the gap in region IV near the IV-V phase boundary 
in the general case, without assuming Zi symmetry. An 
incorrect statement that we made about this in the first 
version of our preprint has been removed. [ E3j also stud- 
ies the shape of the IV-V phase boundary near the SU(4) 
point analytically, showing that it has the opposite cur- 
vature to what is drawn in our Fig. 1. It is possible 
that the phase boundary has zero curvature points close 



to the SU(4) point so that the phase boundary drawn 
in Fig. 1 is roughly correct. On the other hand, there 
are non-zero marginal coupling constants present along 
the phase boundary so that it is difficult to determine it 
accuracy from finite size data and it is possible that the 
boundary drawn in Fig. 1 is quite inaccurate. 

IA would like to thank Rajiv Singh for helpful discus- 
sions. This research was supported in part by NSERC 
of Canada and by the N.S.F. under Grant. No. PHY94- 
07194. 
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FIG. 1. Phase diagram for model Hamiltonian (1.1). In 
Phase I, ground state is composed of fully polarized Ferro- 
magnetic (FM) states for both spin S and orbital T. In phase 
II (phase III), the ground states are antiferromagnetic (AF) 
for 5" and FM fotT (AF for T and FM for S). Phase IV is a 
dimerized phase.u Phase V is a gapless phase. 
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FIG. 2. Renormalization flow for go < and g\ — g-2 < 
according to Eq.(2.26). They flow to zero, the level 1 5(7(4) 
WZW model. We use square root for |go| and \g\\ to show 
the small coupling region clearly. 
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FIG 3. The scaling of the excitation energy 

for 5(7(4) point. The difference between the ground 
state energy and first excited state energy is plotted as 
= (-B(i.i) — -E(o,o))/sin( T ^ T ) vs. 1/lnJ. The fitting line 



(7r/8)(l-l/lnZ). 
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FIG. 4. Energy gap AE and its rate of change dAE(l)/dl 
between triplet-singlet excitation and triplet-triplet excitation 
on x = y line in SU(4) gapless phase. When x = y is much 
smaller than the 5(7(4) point magnitude x — y = 1/4, the 
energy gap AE is bigger. But the rate of change for short 
chain is positive for small deviation from x = y — 1/4 and 
negative for big deviation from x — y = 1/4. This is in 
agreement with the RG flow in Fig. B. 
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Fig. 5 two spin velocities for x=1/4, y=1/8 
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FIG. 5. The scaling of the excitation energies at 

(x,y) = (1/4,1/8) in phase V of Fig. [j] The differ- 
ence between the ground state energy -E(o,o), an d the ex- 
cited state energies E(ij), (S = i, T — j) are plotted as 
Vij = (-B(ij) - £'(o,o))/sin( I ^ r ) vs. 1/lnZ. (v 3 + v )/8vu 
is plotted and reference line at magnitude 0.25 has been 
drawn. Fitting lines shows that v s ~ 0.26, v ~ 0.47, and 
«n = (f s + v )/2 ~ 0.35 in large Z limit. 
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